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ABSTRACT 

Sub-photospheric internal shocks and transverse differences of the bulk Lorentz 
factor in relativistic fireball models of GRB lead to neutron diffusion relative to 
protons, resulting in inelastic nuclear collisions. This results in significant fluxes 
of Vnipu) of ~ 3 GeV and v e {p/) of ~ 2 GeV, scaling with the flow Lorentz factor 
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rj < r/ n ~ 400. This extends significantly the parameter space for which neutrinos from 



inelastic collision are expected, which in the absence of the above effects requires values 



in excess of rj^. A model with sideways diffusion of neutrons from a slower wind into a 
fast jet can lead to production of v^ipp) and v e (p^) in the 2-25 GeV or higher range, 
depending on the value of r/. The emission from either of these mechanisms from GRB 
q at redshifts z ~ 1 may be detectable in suitably densely spaced detectors. 
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1. Introduction 

Neutrinos in three different energy ranges are expected from gamma-ray burst sources. The 
most straightforward and also the hardest to detect are quasi-thermal 10-30 MeV neutrinos 
associated with the stellar collapse or merger event that triggers the burst (Kumar 1999, Ruffert 
& Jahnka 1999). The low energies imply an extremely low detection cross section, which for 
sources at cosmological distances makes the required detector volume prohibitive. At much higher 
energies, e„^2x 10 14 eV neutrinos are expected from the photopion process in jry interactions 
between protons accelerated in shocks and the GRB photons produced by the same shocks 
associated with the fireball. These have a much higher cross section, and significant detection rates 
are predicted (Waxman & Bahcall 1997, Vietri 1998, Rachen & Meszaros 1998). Such neutrino 
events would directly test the shock acceleration mechanism, since they require relativistic protons 
with energies ~ 10 16 eV. A third type of neutrino emission, which is expected in the 5-15 GeV 
energy range (Derishev et al, 1999, Bahcall & Meszaros 2000), occurs only if neutrons are present 
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in the fireball, and the fluences are detectable in km 3 detectors. These neutrinos arise from 
longitudinal decoupling of the n and p flows, and this is expected to occur whether or not shocks 
occur, provided a substantial fraction of neutrons is present in the fireball and the total fireball 
bulk Lorentz factor 77 <; 77^ — 400. 

In this paper we point out that if internal shocks occur below the radiation photosphere, this 
can lead to a regime where rapid diffusion of neutrons both parallel and transverse to the radial 

direction occurs. This leads to neutrinos with energy e u p ~ 3(r//r/ 7r s /j)(2/l + z) GeV for a range 

4/5 

of r\ both below and above a critical value ri^^h ~ w ~ 120. This extends considerably the 
parameter space available for neutrino production from np interactions, in contrast to the smooth 
homogeneous fireball case where n rj n is required. Whereas in typical shocks responsible for MeV 
gamma-rays the neutrino efficiencies from pp are low due to increased coupling and lower optical 
depths (Paczyhski & Xu 1994), here the neutrino fluence is maximized with np collisions in shocks 
at or below the pionosphere, where the inelastic optical depth is <; 1. The flux of such neutrinos 
may be detectable with km 3 Cherenkov detectors currently under planning, for a sufficiently dense 
phototube spacing. Similar or higher neutrino energies and rates are also possible from bursts 
involving a slow wind and a fast jet, where neutrons from the slow wind diffuse sideways into the 
highly relativistic GRB beam. 



2. Sub-photospheric shocks and n diffusion 

In a relativistic GRB outflow consisting of 7, e ,p and n which is initially accelerating with a 
common bulk Lorentz factor T ~ r/r Q , the dynamics is the same for a spherical geometry of for 
an outflow constrained within a channel of angular width 6j <; r -1 . There are two dynamical 
regimes, depending on whether n = L/Mc 2 is greater or less than the critical value 

where £ = n n /n p is the neutron/proton ratio, which for simplicity is taken to be constant over the 
range of interest, ffj~3x 10~ 26 cm 2 is the inelastic cross section near threshold, and r Q = 10 7 r O 7 
cm is the initial radius of the outflow. The parameter rj^ of equation (|l|) can also be written as 
iJtt = f]*( a -K I a T) l ^ A , where ot is the Thomson cross section and 

n* = (7 ^ttTTt) ^ lo3 ^V /4 U + 6" 1/4 (2) 

\4:irr m p c i (l + £) J 

is the critical value of r\ for j,p decoupling in a p, n outflow, which enters also in the definition of 
the photospheric radius (e.g. Meszaros & Rees 2000 for a pure proton fireball). 

For 77 > r/ n , the radial velocity of the n and p becomes decoupled while both components 
are in the accelerating regime where T oc r, leading to relative drift velocities between n and p of 
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Vrei c, and 5-10 GeV neutrinos are expected from inelastic np collisions (Bahcall & Meszaros , 
2000). On the other hand, for 77 < r] n a decoupling of the bulk velocity of the n and p in the 
radial direction is prevented by the fact that the radial motion of both the n and p saturate to 
an asymptotic Lorentz factor F p ~ T n ~ rj at r r s ~ r Q rj before a relative drift v re i — * c has 
been reached, and thereafter the various components coast together. If 77 is constant in time 
and uniform across the channel of angular width 0j, the comoving temperature T' = T {r /r) or 
T' = (T /?])(ry/r) 2 / 3 for r/r Q below or above the value 77, and T Q = 1.2L^t y^ MeV is the initial 
temperature at r Q (henceforth denoting with primes quantities measured in the comoving frame). 
The thermal velocities of the baryons are therefore subrelativistic and the relative CM kinetic 
energy of the n and p is below the threshold energy ~ 140 MeV required for inelastic nuclear 
collisions leading to neutrinos. 

However, for 77 ^ 7/v, it is possible and likely that temporal or angular inhomogeneities in 77 
do exist in realistic GRB outflows. Both types of inhomogeneities imply the presence of density 
gradients, and thus the possibility for neutrons to diffuse into regions whose bulk velocity is 
significantly different from that where they originated. When the neutrons encounter protons or 
other neutrons with relative velocities v re \ <; 0.6c the CM inelastic threshold energy e re i 140 
MeV is reached, leading to pion formation and neutrinos. In order for this to occur with significant 
probability, one requires the inelastic np comoving mean-free path £' np ~ (n^cr,,-) -1 to be comparable 
to the scale of the inhomogeneity in the comoving frame, where n' p = Lj\A. r Kr 2 m v (?Yn[\ + £)] is 
the comoving proton density (henceforth primed quantities are in the comoving frame). In the 
case of internal shocks for which Ar ~ T, which ensures v re i — ► c, the characteristic length scale of 
the shocks is r/T. In the case of sideways diffusion, a thermal velocity vt ~ v re \ — > c also requires 
internal shocks to have re-heated the gas, so both effects are tied to each other, and the transverse 
diffusion is limited to the causal angle T _1 so the lengthscale is again r/T. In the saturated regime 
(which is a prerequisite for internal shocks to occur), when the n,p thermal relative velocities 
v re i — > c so a — > ov, the inelastic optical depth parallel or perpendicular to the shock normal 
becomes smaller than unity above a radius given by 

(rv/r ) = r/ 7r (r/ 7r /? ? ) 3 . (3) 

In the above expression the % dependence includes a factor (1 + £)~ 1//4 , which assumes protons as 
the only targets, for comparison with the radial decoupling case (the case where both n and p are 
targets is given by the same equation (||) by formally setting £ = in the r/^ factor). The internal 
shocks, required to provide v re \ — > c, occur in a coasting regime, at radii 

r s h/r ~ (t v /t )ri 2 , (4) 

where t v is the observer-frame variability timescale on which 77 varies significantly, and 
tv r£ t ~ r /c ~ 10~ 3 s, where t Q is the minimum dynamical timescale in the system over which 
such variation can occur involving a substantial fraction of the total outflow energy. This defines 
a minimum shock radius r s h^ m /r ~ rj 2 , for 7/ < 77*. For radii r s h t m ^ r ^ r w shocks leading to 
Vrei c are possible, while the additional condition > 1 ensures efficient pion production and 
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neutrino emission. This region is shown by horizontal hashed lines in Figure 1. The requirement 
that r s /j m <, r n translates into the condition 

V < rh,.h ~ (to/t v ) 1/5 vt /5 * 120L^ 5 (t o /^) 1/5 , (5) 

For radii above r w , or for values of rj > r) ntS h, shocks and inelastic collisions can also occur in the 
saturated regime at v > T*sh,mi but the mean- free path £ f n is longer than the characteristic system 
length r/r, and the optical depth for inelastic collisions < 1, being equal to the inverse ratio of 
these two lengths. This regime is shown with vertical hashed lines in Figure 1. For rj > 77* > % the 
neutrons have decoupled at a terminal value T n f ~ (3/4)ry 7r (r/ 7r /r7) 1 / 3 (Bahcall & Meszaros 2000, 
Derishev et al 1999), while the protons (which are the component capable of being heated by 
shocks) reach a maximum bulk Lorentz factor T p f ~ 7/* < 7/ (Meszaros &: Rees 2000). Shocks can 
still occur, but the radius at which they occur reaches a maximum value r^^i/To ~ (tv/to)vl — 
constant. 

Below the radius the pionization optical depth t£ <T 1. The neutrino emission from 
inelastic collisions at threshold leading to pion and muon decay contains, on average, one 
Vn and one per neutron, with mean CM energies of 30-50 MeV, and also either a v e or v e 
with an average CM energy of about 30 MeV, resulting from pions and muons decaying at 
rest in the CM (Bahcall & Meszaros , 2000). Typical conditions assumed for internal shocks 
involve collisions between shells with saturated Lorentz factors differing by factors of a few, 
e.g. a root-mean square value rj = ^^a^b = 10 2 r?2 with t]a = 200, t]b = 50 as an example. 
If the luminosity is constant between variations, this corresponds to a mechanical dissipation 
efficiency e = [tja + i]B — ^^VAVb) / {VA + Vb) = 0.2. The typical collision Lorentz factor is 
T rel ~ (1/2)(t/aA/b + Vb/va) - 2.12, and the total CM energy is W = {2m 2 p + 2m p E rel ) 1 / 2 ~ 2.35 
GeV, whereas the CM threshold energy is 2m p + m n ~ 2.01 GeV. Thus, at depths r n > 1 each n 
heated by an individual shock can, on average, collide k w ~ 2 times (and produce a pion each time) 
before its CM energy drops below threshold. The observer-frame energy of the v^v^ (u e D e ) which 
would be predominant in the detector is ~ 50(30)a??/(l + z) MeV, where a ~ 1 near threshold, or 

^ ~ 3a (^ (it;) GeV ■ E - ~ 2 " fc) (ill) GeV • (6 » 

where rj^ S h is given in equation (|5|). The average total v + v CM energy per neutron at 
threshold is e' v9 ~ lOOaA;,,- MeV, which in the observer frame 

GeV. The average neutrino detection cross section averaged over v and v is (Gaisser, 1990) 
a vV ~ 0.5 x 10 _38 (e I/P /GeV)cm 2 , which (per neutron) is a vD ~ 0.6 x W' 37 a(k 7r /2)(i]/ri 7TtSh )(2/l + z) 
cm 2 . The total number of neutrons involved in the shock and diffusion process, with a shock 
dissipation efficiency e = 0.2eo.2, is 

^-(i^)(^)~»^**fcr(i??)- < 7 » 
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Assuming a km 3 detector with Nt ~ 10 S9 NfSQ target protons, and a burst rate out to a 
Hubble radius of 10 3 K h per year, the number of events per year in the detector will be 
Rvv ~ [Nt/ 1 AftD 2 )1Z b N n a ui j, where D is the proper distance out to redshift z. For an Einstein-de 
Sitter universe with Hubble constant H = 65/i65 km/s/Mpc this is 

Ruo,sh ~ 3 E 53 e . 2 a (i^) N t^ n bA (^~~ ~^/TTz) year_1 ' ( 8 ) 

in coincidence with observed GRBs, and independent of the value of rj for rj £j rj^^h, which 
is the limit below which t w <; 1. For values Ti^^h & V the rate (§) must be multiplied 

by ~ (f?7r, s /i/f/) 5 , so between these limits R vi > drops steeply. However, for rj ^> rj n ~ 400, 
even in a laminar flow without any shocks one expects radial decoupling of the np leading to 
inelastic collisions, in which case the rate becomes a factor (&7r£Wi) ~ 2.5 larger, or (Bahcall & 
IVleszaros 2000) Rud ^ 7/year for typical ~ cp„ ~ 10 GeV. 



3. External neutron injection into the beam 

If the progenitor is a collapsing massive star, a relativistic fireball from the collapsed central 
engine, probably in the form of a jet, must punch its way out through the stellar core and the 
envelope in order to become a 'manifest' GRB which emits observable 7-rays. The collision of 
fast baryons from the jet with quasistatic or slow neutrons or baryons from the rest of the star is 
a possible source of high energy neutrinos. We show here that two of the most straightforward 
scenarios have a very low neutrino efficiency, and then go on to discuss a third scenario where the 
efficiency is significantly higher. 

The simplest scenario is that of a beam-dump provided by the outer envelope of a massive 
progenitor star, before the jet has punched through it. The fast neutrons of the jet would 
penetrate ahead of the protons participating in the bow shock, and collide inelastically with 
the upstream baryons at relative Lorentz factors ~ rjj ~ 10 2 — 10 3 . However, a jet with such 
a Lorentz factor takes only t cross ~ r/crj 2 ~ 3 x 10~ 2 ri37/^~ 2 s in the observer frame to reach 
an envelope radius 10 13 ri3 cm, and for a burst of observed duration t& = 30(2/1 + z) s only a 
fraction {t cross /t^) ~ 2 x 10 _3 ri3r/^~ 2 (tft/30)~ 1 of the jet neutrons will collide with envelope targets 
before the latter is punched through. Assuming that the observer is in the jet beam, the isotropic 
equivalent output of jet neutrons is N n = 3 x 10 53 £ , 53?7^ 1 (2^/1 + £). The jet neutrons have a 
source frame energy ej ~ 10 2 ?72 GeV, well above the inelastic threshold for envelope baryons. 
Assuming that the pions decay at rest in the CM, the maximum observer frame energy of the 
Up or Up is e v p ~ 50ary/(l + z) MeV = 2.5ary2(2/l + z) GeV. The average total uv energy 
per neutron is e ~ O.lo;^^/^ + z) GeV and the total detection cross section per neutron is 
o V y ~ 2.5 x 10 -38 aC7r 772(2/1 + z) cm 2 , where £ w is the pion multiplicity. The expected rate of 
events is R uP ,dum P ~ (N t /4TrD 2 )'Jl b a u pN n (t cross /t b ) ~ 1.5 x W~ 2 E 53 a( n r]2 2 r 13 (30s/t b )(2/l + z) 
year -1 , with the same detector and cosmological scalings as in equation @. This is therefore too 
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low for km 3 detectors. 

A second natural scenario involves transverse neutron diffusion into the beam from the 
core of the evolved progenitor star, and illustrates some of the factors affecting neutrino 
efficiency, as well as serving as a preliminary to the third, more efficient scenario discussed in 
the subsequent paragraph. As the progenitor collapse proceeds the core will be heated by a 
shock moving outwards with v = (3 s c to temperatures characteristic of the virial value near 
the central engine region, and it will be permeated by a diffuse bath of thermal neutrinos and 
gamma rays capable of dissociating the metals in the core. The free neutrons then diffuse via 
elastic nuclear collisions, into the jet channel. Taking a core neutron density n n = £ c n c where 
n c = 3 x 10 26 m c r{^ 3 cm~ 3 is the baryon density for a core of m c solar masses, the neutron flux 
diffusing into the channel is cf> n ~ (t>j£/i) 1//2 n n cm~ 2 s -1 after a time t, where £ ~ (n c o'ei)~ 1 ^ an d 
o'el ~ &ir(c/vt) is the elastic cross section, vt ~ (3 s c being the neutron thermal velocity in the 
core. Using an effective transverse area of the channel of A ~ irr 2 6j, the time-integrated fluence 
of neutrons into the channel after time t is N n ± ~ 2tirr 2 6j{3 s (£/l + ^(cfic/a^t) 1 ^ 2 , which is 
N n ,± ~ 10 51 6»_ir 1 1 ^ 2 (m c /2) 1 / 2 (/3 s /3)(2^/l + O(^/30s) 1 / 2 (2/l + z) 1 / 2 , using a core mass m c = 2M , 
an observed duration % = 30(2/1 + z) s and a jet angle 8 = 10 -1 #_i radians. This is comparable 
to the number of baryons going out through the jet, which is N^O 2 ~ 6 x 10 51 E^rj^ 1 2 . i> 
where iVj is the isotropic equivalent output of jet baryons, = 6 x 10 53 £ , 53??^ 1 . The neutrons 
diffusing into the channel collide against the ultrarelativistic baryons in the jet, and as above, 
the maximum observed energy of the or produced is e^ M p M ~ 2.5a?72(2/l + z) GeV 
and the average detection cross section per neutron is a V p ~ 2.5 x 10 -38 0:^^2(2/1 + z) cm 2 . 
However, the transverse inelastic optical depth through the jet encountered by the diffusing 
slow neutrons is T n> ± ~ n n a n r9 ~ 5 x 10 4 L52ry 2 ~ 1 ^io 1 ^-i- Thus only a fraction t~\ of the jet 
baryons can produce 2.5 GeV neutrinos, since the slow neutrons cannot penetrate far into the 
jet. The number of neutrino events in a km 3 detector is then R u p jCor e ~ (Nt /4:irD 2 )7Zbd; / 2Ni } T~ 1 L 
~ 1.5 x 10~ 4 ?72aC7r r io#Ii(2/l + z) 2 year -1 , with detector and cosmological scalings as in equation 
(|8]). This is far too low, and it indicates the critical role played by the transverse optical depth. 

A maximally efficient scenario based on lateral neutron injection can be obtained when the 
jet transverse inelastic depth is of order unity, so that neutrons diffusing inwards can interact with 
the entire fast beam. The transverse inelastic depth of the jet baryons in the lab frame is 

t w± = nja^rO ~ r}\rf l Q(r lr) , (9) 

and unit optical depth is achieved at r/r D ~ 4 x lO^-ir/^i^r^. 1 or about 4 x 10 13 cm for a 
jet with 7/ = 10 3 r73 ~ 1. If the jet is surrounded by a slower outflow, say with T s ~ 10 (and 
T s ^ ®J l )i t ne slow outflow can reach a distance 4 x 10 13 cm in an observer frame burst duration, 
r/cT 2 ~ tb ~ 30 s, and internal shocks in the slow outflow can heat the baryons in it. Inelastic 
collisions will heat the neutrons in the slow flow carried along from lower radii to thermal velocities 
Vt ^ c, and lead also to vv, but their lab-frame energy will be very low, ^ 0.25 GeV. The lab frame 
density in the slow outflow will be 10 2 times larger than in the jet (assuming the same energy 
per solid angle), and neutrons can diffuse into the jet from a depth in the slow outflow which is 
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one-tenth of the width of the jet in the time taken for a neutron to cross the jet, the amount of 
slow mass in this slow width being ten times that in the jet. If the pions and muons are produced 
at rest in the CM, using the CM neutrino energies e' of the previous section the maximum observer 
frame energy of the v or from inelastic collisions with jet baryons is e ~ e'r/(l + -z) -1 



e u ^ ~ 25 ar, 3 J GeV , e UePe ~ 15 a m ^— j GeV , (10) 

where a ~ 1 near threshold. The total observer frame neutrino energy per neutron is 
£u + ~ 0.1qCtt??/(1 + z) GeV, where Cn is the pion multiplicity, and the average detection cross 
section averaged over v and v is a v p ~ 2.5 x lO -3 ^,,- 773 (2/1 + z) cm 2 per neutron. The (isotropic 
equivalent) baryon output in the jet is N b = E/r]m p c 2 = 6 x 10 52 i?53773~ , and the number of 
expected neutrino detection events is 

R»uj ~ (N t /4nD 2 )K b a u 2N b ~ 15 E 53 a^N m Tl b3 hl 5 f - - - _ 




This assumes that the maximum emission occurs at t„.j_ ~ 1, while for t„._l ~ 
i?53r7^~ 1 6'_irj~3 1 6 (tb/30s) _1 (l + z/2) > 1 the rate would be reduced by r~_^. Note that if 
the transverse neutron input N n ± ~ 27r r 2 6 % (3 ' s (ctn s / a ' 2 ~ 1Q 53 ^flsL^rjg-y 1 t^^r^rxa.^-i is 
less than the number of baryons in the jet N b 6 2 = 6 x 10 50 £ , 53^ 2 1^3" , then N n ± must be used in 
equation ([ll|), instead of Nf,; this would be the case only if the wind neutron content is very small, 

The above example used r\ = 10 3 , but the optimum efficiency T n ± ~ 1 can also be obtained 
for other values of rj in such a slow+fast jet model. The rate ([ll]), R V v,j ~ 15/year, is independent 
of the value of 77, since r/ -1 in the number of ejected baryons cancels out with 77 in the neutrino 
detection cross section. The neutrino energies do scale with rj, and the same rate of 15/year 
would be obtained, e.g. for tjj — 10 2 with ^u^u^ ^ 2.5 GeV, 

e v e p e ~ 1-5 GeV, or for r/j — 10 4 with 
e u p ~ 250 GeV, e Ue p e ~ 150 GeV. The latter is an extreme value of r]j, but such energy neutrinos 
would be more easily detectable with first generation km 3 detectors, which could therefore provide 
significant constraints on the jet Lorentz factor. 



4. Discussion 

We have discussed two new mechanisms for producing significant fluences of 2-25 GeV 
neutrinos in GRB fireballs. The first of these involves internal shocks occurring at or below 
the pionosphere (which itself is below the 7-ray photosphere), where the optical depth to 
inelastic nuclear collisions is large. This mechanism requires either temporal or angular 
inhomogeneities (or both) in the relativistic outflow, which otherwise can be considered as 
a simple flow (i.e. characterized by a single mean energy per solid angle and a single mean 
Lorentz factor, aside from order unity variations in it). This mechanism can occur for values 
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of 77 <, % /5 (i /^) 1/5 £ 120L^ 5 r o7 1/5 (t o /t 1 ,) 1 / 5 (2/[l + £]) 1/5 . Thus, it considerably relaxes the rj 
parameter space requirements where neutrinos in the multi-GeV range can be expected. 

The second mechanism involves transverse diffusion of neutrons from a slower outflow into 
a fast jet, at a height where the transverse inelastic optical depth of the latter is close to unity 
for maximal efficiency. This also requires the slow outflow to be heated by internal shocks at 
or below its own pionosphere, and requires a two-zone outflow, characterized by two different 
values of the mean Lorentz factor (or alternatively of the energy per solid angle). This mechanism 
operates optimally at radii where the transverse optical depth is near unity, and predicts nominal 
event rates of ^ 15/year independent of rj and neutrino energies e v & ~ 25 ??3(2/[i + z]) GeV, 
£u e u e ~ 15 ?73(2/[l + z]) GeV which scale with rj, providing a diagnostic for the value of the jet 
Lorentz factor. 

A feature of these two mechanisms is that they can operate even if the outflow has initially 
a low neutron fraction, since inelastic collisions with optical depth ^ 1 tend to convert p into n 
and viceversa with approximately equal probability. They may therefore occur either in neutron 
star mergers or in massive stellar progenitor collapses. For km 3 detectors such as ICECUBE 
(Halzen, 1999), ANTARES (Feinstein 1998), BAIKAL (Beloplatikov et al 1997) or NESTOR 
(Trascatti et al 1998) with sufficiently dense phototube spacing, the predicted event rates are 
Rud ~ 3 — 15/year, which would be detectable above the background in coincidence with GRB 
flashes. However, since these two mechanisms only require sub-photospheric internal shocks (while 
longitudinal decoupling does not require any shocks), 2-25 GeV or higher energy neutrinos could 
be observed from a much wider range of stellar collapse or compact merger events in the Universe 
at z > 1, not necessarily associated with 7-ray flashes. 

This research has been supported by NASA NAG5-2857, the Guggenheim Foundation, the 
Sackler Foundation and the Royal Society. We are grateful to John Bahcall for valuable comments 
and discussions. 
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1 2 3 4 5 lgT] 

Fig. 1. — Radial distance vs. bulk Lorentz factor rj parameter space, showing the pionosphere 
(pionization optical depth 7v = 1), and the photosphere (scattering optical depth r 7 = 1). 
Horizontal hashed regions indicate r n > 1, this being absent above r/ n (where n and p are elastically 
coupled, until decoupling and pionization occur at 7v ^> 1). Vertical hashed regions, only shown 
for part of the range, indicate r n < 1. The bulk Lorentz factor saturation radius r s is common for 
n and p below rj n , while above % protons saturate at r s (p) and neutrons saturate (decouple) at 
r n = 1. Internal shocks are possible above the line r s h(t ), marked by diagonal hashes. 



